Force Modulating Dynamic Disorder: Physical Theory of Catch-slip bond 
Transitions in Receptor-Ligand Forced Dissociation Experiments 
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Recently experiments showed that some adhesive receptor-ligand complexes increase their lifetimes 
when they are stretched by mechanical force, while the force increase beyond some thresholds their 
lifetimes decrease. Several specific chemical kinetic models have been developed to explain the in- 
triguing transitions from the "catch-bonds" to the "slip-bonds" . In this work we suggest that the 
counterintuitive forced dissociation of the complexes is a typical rate process with dynamic disor- 
der. An uniform one-dimension force modulating Agmon-Hopfield model is used to quantitatively 
describe the transitions observed in the single bond P-selctin glycoprotein ligand l(PSGL-l)— P- 
selectin forced dissociation experiments, which were respectively carried out on the constant force 
[Marshall, et at, (2003) Nature 423, 190-193] and the force steady- or jump-ramp [Evans et at, 
(2004) Proc. Natl. Acad. Sci. USA 98, 11281-11286] modes. Our calculation shows that the novel 
catch-slip bond transition arises from a competition of the two components of external applied 
force along the dissociation reaction coordinate and the complex conformational coordinate: the 
former accelerates the dissociation by lowering the height of the energy barrier between the bound 
and free states (slip), while the later stabilizes the complex by dragging the system to the higher 
barrier height (catch). 



Adhesive receptor-ligand complexes with unique ki- 
netic and mechanical properties paly key roles in cell 
aggregation, adhesion and other life's functions in cells. 
A well studied example is the receptors in selectin fam- 
ily which comprises E-, L- and P-selectin interacting 
and forming "bonds" with their glycoprotein ligands. 
These bonds are primarily responsible for the tethering 
and rolling of leukocytes on inflamed endothelium un- 
der shear stress |2() . In particular, in the past two 
years great experimental efforts |3J S Hi LJ 0) have been 
devoted to study the surprising kinetic and mechanical 
behaviors of the bonds between L- and P-selectin and 
P-selectin glycoprotein ligand I (PSGL-1) at the sin- 
gle molecule level: the lifetimes of these bonds first in- 
crease with initial application of small force, which are 
termed "catch" bonds, and subsequently decrease, which 
are termed "slip" bonds when the force increases beyond 
some thresholds. The most important biological mean- 
ing of this discovery is that the catch-slip transitions of 
the PSGL-1— L- and —P-selectin bonds may provide a 
direct experimental evidence at the single-molecule level 
to account for the shear threshold effect in which 

the number of rolling leukocytes first increases and then 
decreases while monotonically increasing shear stress. 

On the theoretical side, it is a challenge to give a rea- 
sonable physical theory or model to explain the counter- 
intuitive bond transitions. Bell IjlOT) firstly suggested that 
the force induced dissociation rate of adhesive receptor- 
ligand complex could be described by, 



fcoff(/) = Cexp[/^ABT], 



(1) 
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where fc^ ff is the intrinsic dissociation rate constant in 
the absence of force, is the distance from the bound 
state to the energy barrier, / is a projection of exter- 
nal applied force along the dissociation coordinate, 
the Boltzmann's constant, and T is absolute tempera- 
ture. The validity of the model has been demonstrated 
in experiments |rH IT^) . Although later at least four 
models have been put forward to explain and under- 
stand various receptor-ligand forced dissociation exper- 
iments fiH IbH n~5^> . they cannot predict catch bonds 
because force in these models only lowers height of the 
energy barrier while shortening lifetimes of the bonds. 
An exception is the Hookean spring model proposed 
by Dembo many years ago \\A Il4f ) . in which a catch 
bond was raised in mathematics. Compared to the ex- 
ponential decay of the lifetimes of slip-bonds with re- 
spect to force in experiments iTL?!) . the model claimed 
that the lifetimes decrease exponentially with the square 
of force In addition, the Hookean spring model 

cannot account for the catch-slip bond transitions in 
self-consistent term. Prompted by the intriguing experi- 
mental observations, three chemical kinetic models have 
been developed. Evans et al. presented a two path- 
ways, two bound states model with rapid equilibrium as- 
sumption between the two states. They suggested that 
the catch-slip bond transitions take place due to applied 
force switching the pathways from the one with slower 
dissociation rate to the fast one (0). Although this in- 
sightful viewpoint well described force jump-ramp exper- 
iments, there are two apparent flaws in physics. First 
if the forced dissociation experiments were performed 
at very low temperatures or higher solvent viscosities, 
force would be independent of the bond dissociations 
since the force in the two pathways and two bound states 
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model only acts on the inner bound states, while these 
states would be "frozen" under this circumstances. The 
other is that force does not accelerate the dissociation 
processes further when the force is sufficiently large for 
the fast dissociation rate is a constant. The next model 
given by Barsegov and Thirumalai with same kinetic 
scheme seems to improve the two flaws in which the dis- 
sociation rates of the two pathways were allowed to be 
force-dependent with the Bell formulas. Unfortunately, 
so many independent reaction constants with arbitrary 
dependencies on the force parameters (total seven pa- 
rameters) and the final dissociation rate depending on 
them in a complicated way make the physical explana- 
tions and the determination of the parameters difficult to 
track. Very recently, a competitive two pathways and one 
bound state model was proposed by Thomas et al. l)l7fl . 
This model is distinct from the others because there is a 
catch pathway therein, which was thought to arise from 
a backward unbinding pathway. But the model is not 
intuitively obvious just like the authors pointed out. 

As one type of noncovalent bonds, interactions of ad- 
hesive receptors and their ligands are weaker. Moreover, 
the interfaces between them have been reported to be 
broad and shallow, such as the crystal structure of the 
PSGL-1— P-selectin bond revealed iftsl) . Therefore it is 
plausible that the energy barriers for the bonds are fluc- 
tuating with time due to either global conformational 
changes or local conformational changes at the inter- 
faces. Association/dissociation reactions with fluctuat- 
ing energy barrier have been deeply studied by statistical 
physicists in terms of rate processes with dynamic disor- 
der lfl9l) during the past two decades. A prototype is the 
ligand rebinding in myoglobin where the rate constant de- 
pends on a protein coordinate lf20|) . Hence it is of interest 
to determine whether the fluctuation of energy barrier 
responses to catch-slip bond transitions. Such studies 
should be meaningful since in the Bell's initial work and 
the other models developed later, the intrinsic rate con- 
stants k® s were deterministic and time- independent, ft 
is possible to derive unexpected results from the relax- 
ation of this restriction. Stimulated by the two consider- 
ations, in the present work we propose that the intrinsic 
dissociation rate in the Bell model is controlled by a con- 
formational coordinate of receptor-ligand complex, while 
the coordinate is fluctuating as a Brownian motion in a 
bound harmonic potential; applied force not only lowers 
the height of the energy barrier as described in Eq.^but 
also modulates the distribution of the conformational co- 
ordinate. In addition to well predicting the experimental 
data, our theory may also provide a new physical mecha- 
nism for the dissociation rates suggested by Bell iflfih and 
Dembo lfl3|) early. 



I. THEORY AND METHODS 

The physical picture of our theory for the forced dis- 
sociation of receptor-ligand bonds is very similar with 



the small ligand binding to heme proteins |20|) : there 
is a energy surface for dissociation which dependents on 
both the reaction coordinate for the dissociation and the 
conformational coordinate x of the complex, while the 
later is perpendicular to the former; for each confor- 
mation x there is a different dissociation rate constant 
which obeys the Bell rate model, while the distribution 
of x could be modulated by the force component along 
x-direction; higher temperature or larger diffusivity (low 
viscosities) allows x variation within the complex to take 
place, which results in a variation of the energy barrier 
of the bond with time. 

There are two types of experimental setups to measure 
forced dissociation of receptor-ligand complexes. First we 
consider constant force mode U 0). A diffusion equa- 
tion in the presence of a coordinate dependent reaction 
is given by 



prcs 



dp(x, t) = D d 2 p 
dt dx 2 
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dx 



p- 



Ox 



ko$(x,f\\)p, (2) 



where p(x, t) is probability density for finding a value 
x at time t, and D is the diffusion constant. The mo- 
tion is under influence of a force modulating potential 
Vf ± (x) = Vi(x) — f±x, where Vi(x) is intrinsic potential 
in the absence of any force, and a coordinate-dependent 
Bell rate. In the present work Eq.^depends on x through 
the intrinsic rate k® s (x), and the distance is assumed 
to be a constant for simplicity. Here f± and f\\ are re- 
spective projections of external force / along the reaction 
and conformational diffusion coordinates: 



fx = /sinfl, 

f\\ = /cos0>O, 



(3) 



and 9 is the angle between / and the reaction coordinate. 
We are not ready to study general potentials here. In- 
stead, we focus on specific V^(a:)s, which make Vf ± (x) to 
be 



V fx (x)=V[x-ri- — )+WU. 



(4) 



where r\ and k are two constants with Length and Force 
dimensions. For example for a harmonic potential 

Vi(x) = V Q + k x (x - x Q ) 2 /2 (5) 

with a spring constant k x in which we are interested, it 
gives 

v ( x - xo -^) =k i( x - x °-it) (6) 



and 



W(fx) = V -f ± x -j±-. 



(7) 



Defining a new coordinate variable y — x — rj — we 
can rewrite Eq. [21 with the specific potentials into 



dp(y,t) 
dt 
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where kf(y) = k s(u + V + f±/ K , f\\)- Compared to the 
original work by Agmon and Hopfield Hfl) , our problem 
for the constant force case is almost same except the reac- 
tion rate now is a function of the force. Hence, all results 
obtained previously could be inherited with minor mod- 
ifications. Considering the requirement of extension of 
Eq. [3 to dynamic force in the following, we present the 
essential definitions and calculations. 
Substituting 



p(y, t) = N exp 



V 



(9) 



into Eq. [SJ one can convert the diffusion-reaction equa- 
tion into Schrodinger-like presentation l)2lf) . 



dt 



D 



.8*4 
dy 2 



U f {y)cj > = -n f ^), 



(10) 



where Nq is the normalization constant of the density 
function at t — 0, and the "effective" potential 



Uf(y) = U(y) + k f {y) 



(11) 



D 



2k B T 



i fav 

2k^T \7)y~ 



d 2 V 
dy 2 



kf(y)- 



We define U (y) for it is independent of the force /. Ea.llUI 
can be solved by eigenvalue technique Ij20j) . At larger D 
in which we are interested here, only the smallest eigen- 
value Ao(/) mainly contributes to the eigenvalue expan- 
sion which is obtained by perturbation approach (|22j) : if 
the eigenfunctions and eigenvalues of the "unperturbed" 
Schrodinger operator 



92 
dy 2 



K = -^ + U(y) 



in the absence of kf(y) have been known, 

n<t>l = -Ktl 



(12) 



(13) 



and kf is adequately small, the first eigenfunction </>o(/) 
and eigenvalue Ao(/) of the operator Tif then are respec- 
tively given by 



Mf) = + #(/) + 



(14) 



E 



U°o(y)kf(y)<f>° m (y)dy AO 



and 



Ao(/) 



A< 0) +A«(/) + A< 2) (/) + •■ 

aSS + I $(y)k f (y)$(y)dy + 

(J <ft(v)kf(vW{v)dy)'' 



(15) 



E 



\0 



AO 



Considering that the system is in equilibrium at the ini- 
tial time, i.e., no reactions at the beginning, the first 
eigenvalue Aq must vanish. On the other hand, because 

$(y)<xexp{-V{v)/2kBT), (16) 

and the square of (f>Q is just the equilibrium Boltzmann 
distribution p eq (y) with the potential V(y), we rewritten 
the first correction of Ao(/) as 



A^(/) = J Pcq (y)k f (y)dy, 
p eq {y) oc exp [-V(y)/2k B T] . 



(17) 



Substituting the above formulaes into Eq. the proba- 
bility density function then is approximated to 



p{y, t) « N exp 



V 



2k B T 



exp[-A o (/)t]0o(/) (18) 



The quantity measured in the constant force experi- 
ments is the mean lifetime of the bond (r) , 



">t Jo 



Q{t)dt, 



(19) 



where the survival probability Q(t) related to the prob- 
ability density function is given by 



Q(t) 



p(x,t)dx = / p{y,t)dy 



exp 



-i(A^(/)+A^ 2) (/))]. (20) 



In addition to the constant force mode, force could 
be time-dependent, e.g., force increasing with a constant 
loading rate in biomembrane force probe (BFP) experi- 
ment (jj) . In principle the scenario would be more compli- 
cated than that for the constant force mode. We assume 
that the force is loaded slowly compared to diffusion- 
reaction process. We then make use an adiabatic approx- 
imation analogous to what is done in quantum mechan- 
ics. The correction of this assumption would be tested by 
the agreement between theoretical calculation and exper- 
imental data. We still use Eq. [21 to describe bond disso- 
ciations with the dynamic force, therefore we obtain the 
almost same Eas. l4lll*l except that the force therein is re- 
placed by a time-dependent function f t . We immediately 
have Q 



cj)(y, t) w exp 



(\ (ft>) + B(t')) dt' 



where the "Berry phase" 



B(t) = I Mft)g- t Mft)dy, 



Mft), (21) 



(22) 



and 4>o(ft) is the first eigenfunction of the time- 
dependent Schodinger operator 



(23) 
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Because the eigenvalues and eigenfunctions of the above 
operator cannot be solved analytically for general kf t , we 
also apply the perturbation approach. Hence, we obtain 
4>o(ft) and Ao(/t) by replacing kf in Eos . 1 1 41 and 1151 with 
k f t . The Berry phase then is approximated to 



B(ft 



E (^r) J<t>°o(y)kf t (y)<t>° m (y)dyx 



}(»)■ 



(24) 



Finally, the survival probability for the dynamic force is 
given by 



Q(t) f» exp 



(\^(f t ,) + \^(f t ,)+B(f t , 



dt' 
(25) 



Different from the constant force mode, data of the dy- 
namic force experiments is typically presented in terms of 
the force histogram, which corresponds to the probability 
density of the dissociation forces p(f) 



P(f) 



dQ 
dt 



'df_ 
dt 



(26) 



Particularly, when the force is a linear function of time 
/ = fo + rt, where r is the loading rate, and zero or 
nonzero of fo respectively corresponds to the steady- 
or jump-ramp force mode in the dynamic force experi- 
ment (|4j), we have 



n/Jo^-^OO + A^Cf) -K/j(.n 



(27) 



exp 



1 



fo 



(a£V) + a£ 2) (/') + £(/'))<#' 



II. RESULTS 



We consider a bounded diffusion in the harmonic po- 
tential Eq. [5] Then H reduces to a harmonic oscillator 
operator with 



*T, -, Dk x ( k x y 2 
U(y) = I - 1 



(28) 



2k B T \2k B T 
Its eigenvalues and eigenfunctions are 

A° = nDk x /k B T (29) 

and 

4P n {z) = 2-/ 2 7 r- 1 / 4 (n!)- 1 /2 e - 2 /2 jffn(z); (30) 

respectively, where z = (k x /2k B TY^ 2 y and H n (z) is the 
Hermite polynormials 1)221) . Given that the intrinsic dis- 
sociation rate satisfies the Arrcnhenius form 



k° oS (x) = k exp[-AG*{x)/k B T] , 



(31) 



where the height of the energy barrier along the reaction 
coordinate AG'(i) is a function of the conformational 
coordinate x. According to the form of barrier, we first 
analyze two simple and meaningful cases. 
Bell-like forced dissociations. The simplest function 
of the energy barrier might be linear with respect to x, 



AG* (a;) = AG;?, + k g (X — Xo), 



(32) 



where AGq is the height at position xo, and the slope 
kg > (its dimension Force) for the perturbation re- 
quirement in solving Eg. 1101 According to Eas. 1151 and 
1241 we easily get 



K (/) = fcoexp 



-/3AG, 



exp 



-k 2 
(3Dk x 



■ exp 



exp 



2/3AGJ 
-> i ( i ./ - ^/j 



(33) 



^— ' r).r),l \ 



and 
B(/t) 



^/t|| 



(3D 2 k 2 



exp 



-2/3AG* + 



(34) 



exp 



2/3 ev. 



Ml 



g f 
-Jt± 



oo 



where /3 = l/k B T. For large 13 or fc 2 (or very small T), 
the second correctness and the Berry phase tend to zero. 
Under these limitations the first eigenvalue of Eq. [21 is 
approximated to be 



Ao(/) ~ fcoexp 



(3k 2 a 



Here we define a new distance 



d t =£* cos 6< -(sin( 



exp[/9d*/]. (35) 



(36) 



where £ = k g /k x whose dimension is Distance. We see 
that the presence of the complex conformational coordi- 
nate could modify the original Bell model in novel ways: 
(i) cF > 0, Eq. |23 is indistinguishable from the origin 
Bell model, although the projection distance cos 8 from 
the bound state to the energy barrier may be increased 
or decreased in terms of the orientation of the applied 
force. In particular, if the force is antiparallel to x, i.e., 
sin = — 1, we get a Bell- like rate expression with a "dis- 
tance" C; (ii) ds — 0, the force does not affect disso- 
ciations of the bonds, which have been named "ideal" 
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bonds by Dembo (d Q); (iii) £ < 0, the force slows 
down dissociations of the bonds. It is "catch" bonds in 
which we are interested. In contrast to the catch behavior 
suggested by Dembo (0), the rate decays exponentially 
with respect to the force instead of the square of the 
force. Given the linear function Eo . 1321 and f± > 0, in- 
creasing of the force only stabilizes the bonds by dragging 
the system to the higher energy barriers (catch) , whereas 
the other force component f\\ destabilizes the complex by 
lowering the energy barriers (slip). Therefore the sign of 
the distance cfi in fact reflects a competition of the two 
contrast effects of the same force. 

Dembo-like forced dissociations. Another function 
of the energy barrier is a harmonic with a spring con- 
stant k„ 



AG*(x) = AG* + k g (x - x{) 2 /2, 



(37) 



Because 

(2) 



where AG^ is the barrier height at position x\ 

for any form of the barrier height, the dependence of A 
and B(f t ) on D is the same from Eq.|551 we only consider 
the large D limitation in the following. Hence we have 



Ao(/) « k 



k x + kq 



exp 



f3k g (f± - k x (xi - x )) 2 



2k x (k x + kg) 



(38) 



Given sin# > and x\ — xq > 0, we find that there is a 
interesting transition from slip to catch bond when the 
force increase over a threshold k x (x\ — xq)/ sin#; other- 
wise only catch bond presents. We note that the latter is 
very similar to the result proposed by Dembo |)13() even 
their physical origins are completely different: both of 
them exponentially dependent on the square of the force. 
Comparison with the experiments. In the constant 
force rupture experiment of the PSGL-1— P-selcctin com- 
plex, the dissociation rate as the inverse mean lifetime of 
the complex first decreased and then increased when the 
applied force increased beyond a force threshold (0). We 
now can easily understand this counterintuitive transi- 
tion according to the previous discussion about the Bell- 
like dissociation rate: the dissociation effect of /|| regains 
its dominance when the force is beyond the threshold. 
Although in principle we can construct various barrier 
height functions which result into catch-slip transitions, 
the most simplest form may be a composition of two lin- 
ear functions 



AG* (a 



= f AGl(x) = AG* 
(AG*(x) = AG* 



k c {x 
k s (x 



Xb), x < x b 
Xb), X > x b 



(39) 



where we require that the distances defined in Ea. l36l with 
k c and k s are respectively minus and positive. For con- 
venience, their absolute values are correspondingly de- 
noted by d\ and d\. Define two "intrinsic" dissociation 
constants 



k exp[ 
-- fc exp[ 



-PAGl(x )], 
-0AQl(x o )]. 



(40) 



AG, 



ag| 



AG 



| AG* 


/ 

/ Jill ^eq 


fx 






> 








> 



X 







FIG. 1 Schematic diagram of the height function of the energy 
barrier with respect to the coordinate x. AGj and AG| are 
the value of the linear functions in Ea. Kffll at position xo, while 
AGj is the intersection of the functions at Xb- The bold solid 
line is the minimum model which is used to fit the experiment. 
The shaded area represents the equilibrium distribution of the 
conformational coordinate under the potential V(y) ( Eg, 1171 . 



Fig- CI shows the characteristics of the function. We then 
have 



Ao(/) 



exp 



2 

xerfc 
k, 



exp [-0 c f] 



A &C 

A + v 



(ik x (3 

h f \ sin ( 

2 J V 2k* 



fexp 



exp [pSj] 



(41) 



xerfc 



A + 



(3k x 
2 



where A = Xb—xo, and the complementary error function 



erfc(a;) 



" x dx. 



(42) 



Before fixing numerical values of the parameters in 
Eq. 1411 we first simply analyze the main properties of 
Ao(/) given A ^> 0: (i) in the absence of force, due to 
erfc(— oo) = 2 and erfc(+oo) = 0, we have 



X ^k c eMPk 2 c /2k x ), 



(43) 



which is the same with that obtained by Agmon and 
Hopfield (ii) if force is nonzero and smaller, 



A « k c o exp{/3k 2 /2k x )exp(-0 c f), 



(44) 



which means that the bond is catch; and finally (iii), 
when the force is sufficiently large, Eq.^] reduces to 



Ac w k s exp(pk 2 /2k x )exp(pdlf). 



(45) 
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TABLE I Comparison of the parameters of the present the- 
ory, and the two-pathw ay a nd one energy-well model pre- 
sented by Thomas et al. (u7l) on the constant force (cf) and 
the dynamic force (df) modes. The parameters for the slip 
behavior of the P-selectin are also listed as a reference (123 





d c nm 


d s nm ) 


s§ sec 


&o sec 1 


Experiment 




0.14 




0.2 


Dynamic disorder by us 


1.2 


0.22 


23.2 


1.68 


Two-pathway one-well (cf) 


2.2 


0.5 


120 


0.25 


Two-pathway one-well (df) 


0.4 


0.2 


20 


0.34 



It is the ordinary slip bond. 

There are total eight independent parameters present- 
ing in Eq. ED 6, k x , k c , k s , A, £*, k , and AGl(x ). 
It is not necessary to determine all of them, which is 
also impossible only through fitting to the experimental 
data ||3J). For example, the latter two parameters are 
lumped into k$, while and always presents with cos0 
together. What we really concern with is the coefficients 
of the force and the factors before the error functions 
in Eq. They can be obtained by least square fit. 
Guided by the properties Easl4^1andl4*4l the fitting pro- 
cess in fact is simple. Even so, we are still able to fix 
all parameters from the fitting results if we study a min- 
imum model in which the slop k s is zero and 9 = tt/6. 
Here the particular value of the angle is actually of no 
particular significance and it is only as a reference. We 
immediately have: k c ss 0.60 pN, k x w 0.21 pN nm -1 , 

~ 0.25 nm, A w 33 nm; the other interesting parame- 
ters see Tab. 1, where the values are independent of the 
angle. Substituting these values into Eg. 1411 and accord- 
ing to Eq.^3 we calculate the mean lifetime of the PSGL- 
1— P-selectin complex with respect to different constant 
force in Fig. [5] the agreement between theory and the 
experimental data is quite good. 

More challenging experiments to our theory are the 
force steady- and jump-ramp modes 10). Under large D 
limitation, Eq. [23 reduces to 



P(fJo) 



Mf) 



■ exp 



1 



Ao (/')#' 



So 



(46) 



We see that the mean lifetime can be extracted from 
the above equation by setting / = /q, i.e., (r) = 
l/rP(/o,/o)- We calculate the dissociation force distri- 
butions of the steady- and jump-ramp modes at three 
loading rates to compare with the BFP experiments per- 
formed by Evans et al. (4). Here we are not ready to 
fit the experiments afresh; instead we directly apply the 
parameters obtained from the constant force mode to 
current case. Because the BFP experimental data is for 
dimeric ligand PSGL-1, whereas our parameters are from 
monomeric ligand sPSGL-1. Therefore it is necessary 
to map our predictions for the single bond sPSGL-1— P- 




20 30 

Force (pN) 

FIG. 2 The mean lifetime as a function of force for the bonds 
of dimeric P-selectin with monomeric sPSGL-1 (square sym- 
bols) Q and the rescaled dimeric PSGL-1 (circle symbols) 
from Ref. ljl7j) . The two dash curves are respectively cal- 
culated by the two addition terms in Eq. 1411 with the same 
parameters. 



sclcctin to the double bonds. A natural assumption is 
that the two bonds share the same force and fail ran- 
domly. The same assumption has been used in previous 
works ifS ITtI) . Hence, the probability density of the dis- 
sociation force for the double bond PSGL-1— P-selectin 
complex is related to the single case by 



P d (f,f )=P(f/2,f /2Y 



(47) 



Fig. presents the final result. We see that the theoret- 
ical prediction agrees to the data very well. Hence we 
conclude that the adiabatic approximation proposed at 
the beginning is reasonable. The previous works 10 Il7|) 
have claimed that they could not fit the experimental 
data from the constant force experiment using atomic 
force microscopy (AFM) and the force jump-ramp ex- 
periment using BFP with the same parameters, e.g., see 
Tab. 1. The authors simply contributed it to the dif- 
ferent equipment and biological constructs though the 
experiments studied the same complexes Q 0) . Our cal- 
culations however show that the mechanical parameters 
defined by us have almost the same values. In addition, 
the tendencies of our density functions for the first two 
panels of the second array in Fig. are closer to the data 
than that predicted by the two-pathways models if^ llTl) . 

The density functions and the force histograms in the 
experiments reach the maximum and minimum at two 
distinct forces, which are named / m i n and / ma x in the 
following, respectively. This observation could be under- 
stood by setting the derivative of Eq. E3with respect to 
/ equal to zero, 



(48) 
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FIG. 3 The probability density of the dissociation forces 
PdCfi/o) under the different loading rates predicted by our 
theory (solid curves) for the PSGL-1— P-selectin complex. 
The symbols are from the force steady- and jump- ramp exper- 
imental data (4). The apparent deviations between the theory 
and the data in the last column may be from the invalidation 
of the assumption of two equivalent bonds at higher loading 
rates. 



We immediately see that the values of / m i n and / max 
must be larger than the catch-slip transition force f c 
for the left term in Eq. 0H1 is negative as the bond is 
catch. Indeed the experimental observations show that 
the force values at the minimum histograms are around 
a certain values even the loading rates change 10-fold, 
(see Figs. 2 and 4 in Ref. The above equation could 
have no solutions when the loading rate is smaller than 
a critical rate r c , which can be obtained by simultane- 
ously solving Ea. 1481 and its first derivative. We estimate 
r c w 6 pN/s using the current parameters, while the force 
/max = /min ~ 13 pN (about 26 pN in the double bond 
cases). If r < r c , then the density function is monotonous 
and decreasing function. Therefore, the most probable 
force at the bond dissociation is zero. The most inter- 
esting characteristics of Eq.2Hl are the dependence of the 
maximum and minimum forces on the loading rate. In 
particular the latter is an important index in dynamic 
force spectroscopy (DFS) theory since it corresponds to 
the most possible dissociation force |24j). When the load- 
ing rate is sufficiently large, and correspondingly / max is 
larger, the approximation of Eq.^]at large force Ea.1441 
implies that 



/max ~ f-j j In 



1 

(3d s 



(3rd\ 



fcgexp \J3k 2 s /2k x 



oc lnr. 



(49) 



The experimental measurement supported this predic- 
tion; see Fig. 3A in Ref. (Q). On the other hand, due to 
that / m in is very close to f c at the larger loading rate, 



employing the Taylor's expanding approach we have 

Ao(/) = AoC/c) + l^(fc)(Sf) 2 + o[(Sf)% (50) 
where Sf = f — f c . Substituting it into Eq. 0SJ we get 



/min ~ fc + Aq(/c) J r ~^-{fc) <* 



(51) 



It means that / m ; n tends to f c very fast. Different from 
/max, the loading rate dependence of / m in is an intrinsic 
property of the catch-slip bond; a unique requirement 
is that the dissociation rate Ao(/) has a minimum at 
the transition force f c . Therefore / m inS observed in ex- 
periment performed by Evans et al. (jjj) are almost the 
catch-slip transition force observed in the constant force 
rupture experiment performed by Marshall et al. (0). In- 
deed, the force values of the minimum force histograms 
for the former are about 26 pN, while the transition force 
for the latter (dimeric PSGL-1— P-selectin) is also about 
26 pN. We know that the dissociation forces distribution 
of a simple slip bond only has a maximum at a certain 
force value that depends on loading rate |25Tl. There- 
for the catch-slip bond can easily be distinguished from 
the slip case by the presence of a minimum on the den- 
sity function of the dissociation forces at a nonvanished 
force. Because the above analysis is independent of the 
initial force /o, in order to track the catch behaviors in 
the force jump-ramp experiments, /o should be chosen to 
be smaller than / c . 



III. CONCLUSIONS AND DISCUSSION 

Compared to the chemical kinetic schemes, our the- 
ory should be more attractive on the following aspects. 
First of all, we suggest that the counterintuitive catch- 
slip transition is a typical example of the rate processes 
with dynamic disorder. Because this concept has been 
broadly and deeply studied from theory and experiment 
during the past two decades, extensive experience and 
knowledge could be used for reference. For example, we 
suggest that a new receptor-ligand forced dissociation ex- 
periment could be performed over a large range of tem- 
peratures and solvent viscosities. According to Eq. if 
the viscosity is so higher that D — > 0, we could predict 



p(x,t) mp(x,0) exp[-ifc off (a;, /||)]. 



(52) 



We know that such a dissociation reaction is a typical ex- 
ample of the rate processes with static disorder <[li^. In 
addition that the survival probability of the bond con- 
verts into multiple exponential decay at a single force 
from the single exponential decay at the large D limita- 
tion (see Eq. 120(1 , the mean lifetime is 



p(x, 0)kj(x, /||), 



(53) 
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which means that the catch-slip bond changes into slip 
bond only. Then our theory gives a intuitively obvious 
physical explanation of catch bonds in an apparent ex- 
pression fEa. I35f) : they could arise from a competition of 
the two components of applied external force along the 
dissociation reaction coordinate and the molecular con- 
formational coordinate; the former accelerates the dis- 
sociation by lowering the height of the energy barrier, 
while the latter stabilizes the complex by dragging the 
system to the higher barrier height. Finally, the time- 
dependence of the forced dissociation rates could be in- 
duced by either global conformational changes of the 
complex or local conformational changes at the interface 
between the receptor and ligand; no separated bound 
states and pathways are needed in the current theory. 
Therefore it is possible that one cannot find new sta- 
ble complex structures through experiments or detailed 
molecular dynamics (MD) simulations. 

Even there are many advantages in the present theory. 
We cannot definitely distinguish which theory or model 
is the most reasonable and more close real situations 
with existing experimental data. Moreover, except the 
coarse-grain physical picture our theory does not reveal 
the detailed structural information of the catch behavior 
of the ligand-receptor complexes, while biologists might 
be more interested in it. We could correspond the in- 
creasing height of the energy barrier with respect to the 
conformational coordinate to the hook structure l|2q) or 
more affinity bound states l)27j) , however we believe that 
further single-molecule experiments including microma- 
nipulation experiments and fluorescence spectroscopy, 
more crystal structure data and detailed MD simulations 
from the atomic interactions are essential to elucidate 
the real molecular mechanism of the catch bonds. 



Fei Ye for 
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